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Abstract — In recent years, some spectrum sensing algorithms 
using multiple antennas, such as the eigenvalue based detection 
(EBD), have attracted a lot of attention. In this paper, we are 
interested in deriving the asymptotic distributions of the test 
statistics of the EBD algorithms. Two EBD algorithms using 
sample covariance matrices are considered: maximum eigenvalue 
detection (MED) and condition number detection (CND). The 
earlier studies usually assume that the number of antennas (K) 
and the number of samples (JV) are both large, thus random 
matrix theory (RMT) can be used to derive the asymptotic 
distributions of the maximum and minimum eigenvalues of the 
sample covariance matrices. While assuming the number of 
antennas being large simplifies the derivations, in practice, the 
number of antennas equipped at a single secondary user is usually 
small, say 2 or 3, and once designed, this antenna number is fixed. 
Thus in this paper, our objective is to derive the asymptotic 
distributions of the eigenvalues and condition numbers of the 
sample covariance matrices for any fixed K but large N, from 
which the probability of detection and probability of false alarm 
can be obtained. The proposed methodology can also be used 
to analyze the performance of other EBD algorithms. Finally, 
computer simulations are presented to validate the accuracy of 
the derived results. 

Index Terms — Spectrum Sensing, Cognitive Radio, Random 
Matrix Theory. 

I. Introduction 

Spectrum sensing is one of the key elements in opportunistic 
spectrum access design |fl~), El- In literature, there are a 
few main categories of practical specnum sensing schemes 
for OSA including^: energy detection, cyclostationarity based 
detection, and eigenvalue-based detection (EBD); see, e.g. 
O and references therein. Energy detection is simple for 
implementation, however, it requires accurate noise power 
information, and a small error in that estimation may cause 
SNR wall and high probability of false alarm (4). For cyclosta- 
tionarity based detection, the cyclic frequency of PU's signal 
needs to be acquired a priori. On the other hand, EBD 0, 
first proposed and submitted to IEEE 802.22 working group 
0, performs signal detection by estimating the signal and 
noise powers simultaneously, thus it does not require the cyclic 
knowledge of the PU and is robust to noise power uncertainty. 

In general, two steps are needed to complete a sensing 
design: (1) to design the test statistics; and (2) to derive 

'Matched-filtering method requires perfect knowledge of the PU's signal 
received at SU, which is almost impossible in cognitive radio scenario. Thus 
matched-filtering method is not treated a practical spectrum sensing scheme 
here. 



the probability density function (PDF) of the test statistics. 
For spectrum sensing with multiple antennas, the test statis- 
tics can be designed based on the standard principles such 
as generalized likelihood ratio testing (GLRT) (7j, or other 
considerations 0, 0, 0, iflOl . Surprisingly these studies 
all give the test statistics using the eigenvalues of the sample 
covariance matrix. It is thus important to derive the PDF of the 
eigenvalues so that the sensing performance can be quantified. 

For EBD schemes, the PDF of the test statistics is usually 
derived using random matrix theory (RMT); see, e.g., 0, 0, 
0, IflOl . In fact, the maximum and minimum eigenvalues of 
the sample covariance matrix have simple explicit expressions 
when both antenna number (K) and ample size (N) are 
large ifTTI . iFFD . It is reasonable to assume that TV is large 
especially when the secondary user is required to sense a weak 
primary signal, in practice, however, the number of antennas 
equipped at a single secondary user is usually small, say 2 
or 3. Thus the results obtained under the assumption that 
both K and N are large may not be accurate for practical 
multi-antenna cognitive radios. In this paper, our objective is 
to derive the asymptotic distributions of the eigenvalues of the 
sample covariance matrices for arbitrary K but large N. The 
asymptotic results obtained form the basis for quantifying the 
PDF of the test statistics for EBD algorithms. 

It is noticed that there are studies on the exact distribution 
of the condition number of sample covariance matrices for 
arbitrary K and N lfl3l . the formulas derived however are 
complex and cannot be conveniently used to analyze the sens- 
ing performance. Furthermore, there are no results published 
for the case when the primary signals exist. 

The rest of the paper is organized as follows. Section Gj] 
presents the system model for spectrum sensing using multiple 
antennas. Two EBD algorithms are reviewed in Section [TTTJ 
including maximum eigenvalue detection (MED) and condi- 
tion number detection (CND) algorithms. In Section IIV1 we 
derive the asymptotic distributions of the test statistics of the 
two EBD algorithms for the scenario when the primary users 
are inactive. In Section [V] the results are derived for the 
scenario when there are active primary users in the sensed 
band. Performance evaluations are given in Section [VI] and 
finally, conclusions are drawn in Section IVIII 

The following notations are used in this paper. Matrices 
and vectors are typefaced using bold uppercase and lowercase 
letters, respectively. Conjugate transpose of matrix A is de- 



noted as A H . ||a|| stands for the norm of vector a. E[-] denotes 
expectation operation; — ^— ^ and = stand for "convergence 
in distribution", "almost surely convergence"; and "defined 
as", respectively. 

II. System Model 

Suppose the SU is equipped with K antennas, which are 
all used to sense a radio spectrum for collecting N samples 
each. One of the following two scenarios happens. 

Scenario 1 (Si): There are t (t > 1) active PU trans- 
missions in the sensed band, and the sampled outputs, x(n), 
n = 0, • ■ ■ , N — 1, can be represented as 

x(n) = Hs(n) + u(n), (1) 

where u(n) is the (A" x 1) noise vector, s(n) is the (t x 1) 
signal vector containing the transmitted signals from the active 
PUs, and H = [h l7 • ■ • , h t ] is the (Kxt) channel matrix from 
the active PUs to the SU. 

Scenario (So): There are no active PUs in the sensed 
band, thus the sampled outputs collected at SU are given by 

x(n) = u(n). (2) 

We make the following assumptions for the above models: 

(i) The noises u(n) are independent and identically dis- 
tributed (iid) both spatially and temporally. Each element 
follows Gaussian distribution with mean zero and vari- 
ance (J 2 ,. 

(ii) For each n, the elements of s(n) are iid, and follow Gaus- 
sian distribution with mean zero and variance a 2 s . Thus 
the covariance matrix of s(n) is R s = ¥,[s(n)s H (n)] = 
all. 

(iii) s(n) are u(n) are independent of each other. 

(iv) The eigenvalues of HH^ are not all identical. 

Two cases of Gaussian distribution are considered: (1) Real- 
valued case: we denote u(n) ~ Af(0, all) and s(n) ~ 
J\f(Q, a^T); (2) Complex-valued case: We denote u(n) ~ 
CAf(0,alI) and s(n) ~ CAf(0,a 2 s l); 

When there are active PUs, we define the average received 
signal-to-noise ratio (SNR) of PUs' signals measured at SU 

spt ilk 112^2 

as : SNR = % £ J ° ■ Given the received samples, X = 
[x(0), x(l), • • • , x(7V — 1)], the objective of spectrum sensing 
design is to choose one of the two hypothesis: Hq: there are 
no active primary users in the sensed band, and r H.\. there 
exist active primary users. For that, we need to design the test 
static, T(X), and a threshold e, and infer %i if T(X) > e, 
and infer H if T(X) < e. 

III. Eigenvalue based Detections 

Let us define the sample covariance matrix and covariance 
matrix of the measurements x(n) of the sensed band as 

R x t -XX fl = -^x(,)x ff ( n ), (3) 

n=0 

K x = E[x(n)x H (n)], (4) 



respectively. For a fixed K, under Scenario So and when N — » 

00, we have 

R x -»■ R x = <t 2 I, (5) 

which means that all the eigenvalues of the covariance matrix 
are equal to a\. However, under S\ and when N — s- 00, R x 
approaches 

R.t = a 2 s lill H + a 2 u l. (6) 

Based on assumption (iv), the eigenvalues of © can be 
ordered as p\ > P2 > ■ ■ ■ > pk > &l- 

Denote the ordered eigenvalues of the sample covariance 
matrix R x as Ai > A2 > • • • > \k- We consider the following 
two EBD algorithms. 

(1) Maximum Eigenvalue Detection (MED): For MED, the 
test static is chosen as J8): 

T M(X) = ^. (7) 

(2) Condition Number Detection (CND): The CND 

chooses the following test statistic J5]: 

T (c) (X) = b-. (8) 

An essential task to complete the sensing design is to 
determine the test threshold, which affects both probability 
of detection and probability of false alarm. To do so, it is 
important to derive the PDF of the test statistics. 

For arbitrary (K, N) pair, the closed-form expressions of 
the PDF of the test statics are in general complex |[T3l . 
In practice, the primary users need to be detectable in low 
SNR environment. For example, in IEEE 802.22, the TV 
signal needs to be detected at -20dB SNR with 90% target 
probability of detection and 10% target probability of false 
alarm. To achieve that, the number of samples, N, required 
for spectrum sensing is usually very large. In the paper, we 
thus turn our attention to derive the PDFs of the test statistics 
for any fixed K, but large N. The asymptotic distributions 
(when N — s- 00) of the test statics of the EBD algorithms will 
be derived for Sq and Si scenarios, respectively. 

IV. Asymptotic Distributions of the Test 
Statistics under Scenario S 

Let Ai = £t for i = 1 , • • • , K, and denote A = 

diag{Ai,A 2) --- "a^}. Define ft = Va(A 4 - 1) for i = 

1, -- - ,K, and B = diag {/3 X , fc, ■ ■ ■ ,f} K }. Note that ft > 

P2>--->PK- 

Under So, by Theorem 1 and (2.12) in lfl4ll . we have the 
following proposition. 

Proposition 1: For real-valued case, when N — > 00, the 
limiting distribution of B, giciPx , ft , • • • , ftf), is given by 

= ^(A^xpU-^ft 2 ) J] (ft -ft), (9) 

V i=l / l<i<j<K 



where 



Ci(K) 



2^+ 3 )/ 4 nf = ir[i(A' + i- l )] 



(10) 



Here T(z) = J °° t z x e t dt is the Gamma function. 
For real-valued case, and when K = 1, we have /3i 

W(0, 2), i.e., gi {p x ) = ^e-T. When K = 2, we have 



flaC8i,j8a) 



.£i±£ 



2 5 /2, 



{p\-fa)e 4 , (11) 



For complex-valued case, by Lindberg's central limit theo- 
rem 1151 . Vn(^t — I) converges to a Hermitian matrix with 
elements above diagonal being complex Gaussian distribution, 
a so-called Gaussian unitary matrix. Thus, according to the 
joint density of the eigenvalues of a Gaussian unitary matrix, 
similar to Theorem 1 in 1 14 1 we have the following: 

Proposition 2: For complex-valued case, when N — > oo, 
the limiting distributions of B, gK(Pi,(32, ■ ■ ■ ,Pk), is given 
by 



9k{Px,--- ,0k) 
C 2 {K) cxp 



\hf) n ( 

i=\ J l<i<j<K 



(A-ft-) 2 ,(i2) 



where 



C 2 (K) = Kl{2n) 



-K/2 



8 r (i + i)' 



(13) 



For complex-valued case, when K = 1, we have g((3i) 
7=e - ^, i.e., /3i - W(0, 1). When K = 2, we have 



ff 2 (/3i,/32) = ^(/3i-/3 2 ) 2 e- 



(14) 



For a given K, the limiting distribution, g^^x), of fa, i = 
1, • • • , K, can be calculated from (0 or 1121 . 

A. MED 

For MED, T^(X) = |j- = Ai, thus \ZF(T( m )(X)-l) = 
/3i. We have the following: 

Theorem 1: Under So, for a fixed K and when iV — > oo, 



(r (m) (x) - 1 



J (0) 

t _fS" ! 



(15) 



whose distribution, fx (%), is the same as the limiting 
distribution of <7k,i(#). 

B. CND 

For CND, T( C )(X) = h- = The following theorem 
states the limiting distribution of the condition number of the 
sample covariance matrix. 

Theorem 2: Under So, for a fixed K and when N — > oo, 



TV (V c >(x) - l) 



d (0) 



-if ' 



(16) 



whose PDF and cumulative density function (CDF) are given 
by 



F ( k\x) = 



9k{$k + x,(3 K )df3 K , 

Pk+x 



(17) 

g K (Pi,P K WidpK, (18) 

-oo J p K 

respectively, where is the joint distribution of (3\ 

and /3k- 

Proof: See US. ■ 
Next, we derive the closed-form expressions of the limiting 

distributions of the condition number of the sample covariance 

matrix for K = 2 case. 

Lemma 1: Consider Scenario So- For K = 2 and N —> oo, 



y/N (V (c) (X) - 1 



D 



whose PDF is 



#(a;) = i a; exp(-x 2 /8), 



for real-valued case, and 



JO) 



x > 0, 



x > 0. 



(19) 



(20) 



(21) 



20? 

for complex-valued case. 

C. Threshold Determination 

From the above two subsections, it is seen that for a 
given K and when N — > oo, the regulated test statistics, 



JV(Jw(X) — 1) converges to a random variable with PDF 
fx(x) and CDF F^\x), where i = m for MED and i = c 
for CND. Here, we look at how to set the decision threshold 
e to achieve a target probability of false alarm, Pf. Since 



P f = Prob (T«(X) >e|V 

= Prob f\/iV (t (?) (X) - l) > V^V(e 

/]V(e-l) 



0l<5o) 



= 1 - 



/£>(t)dt 



l-F®(VN(e-l)). 



Thus the decision threshold is determined by 



1 



[F^]-\l-P f 



(22) 



(23) 



where x = [F^-] denotes the inverse function of y = 



V. Asymptotic Distributions of the Test Statistics 
under Scenario Si 

In this section, we derive the asymptotic distributions of 
the test statistics under Si. Consider the ordered eigenvalues 
of covariance matrix R^, p\ > p 2 > ■ ■ ■ > pk, and the 
ordered eigenvalues of sample covariance matrix R T , Ai > 
Aa > • • ■ > Ax . With the same notation as in fl4l . let the 



multiplicities of the eigenvalues of H x be qi,qz,- 
is 



, q r . That 



pi = ■ ■ ■ ■■ 

Pqi+l = ■ 
PK-q r +l 



Pqi = Ml) 

' = Pqi+q 2 -l = M2, 
--■■■= PK = Mr, 



where 



Moreover, let 



Mi > M2 > • • • > Mr 



r ff R,r = a, 
E ff r ff R 2 rE = b, 



where 



A 

B = 



= diag(pi,p 2 , 
diag(Ai, A 2 , 



,Pk), 



and r and E are corresponding eigenvector matrices. We then 
partition the matrices A and B into block matrices as follows: 







• 


• °\ 







M2I • 





A = 










v 


• 


• Mrl/ 




fB 1 


• • 


• °\ 







B 2 • 


■ 


B = 










• ■ 





(25) 



Note in both A and B, the kth diagonal sub-block is qj. x q^ 
matrix. 

Following the proofs of Proposition 1 and Proposition 2, 
and Theorem 1 of lfT4l . we have the following results. 

Theorem 3: The limiting distribution of \fN (j^- — I J is 
given by (0 and $Y2\ for real-valued case and complex-valued 
case, respectively, with the parameter K being replaced by g^. 
Furthermore, y/W ( — — I ) is asymptotically independent of 



Let 71 = VW(^ - 1) and lK = VN{^- - 1). From 
Theorem[3] 71 and jk are asymptotically independent, and the 
limiting distributions of 71 and 7^ are g qi ^{x) and g qi ,,q T {x), 
respectively, which are defined in Section [IV] 



A. MED 

For MED, T( m )(X) 



_ Ai 



, thus ^ = a( Tn) T('")(X) where 



a 



(m) 



We have the following theorem. 



Theorem 4: Under <Si, when N - 



■ 00, 



'if ' 



(26) 



whose PDF is (7 gi! i(x). 

If there are no repeated maximum eigenvalues (i.e., qi = 1), 



A/"(0, 2) for real-valued case, and m^ 5 - 7^(0, 1) for 



complex-valued case. 




Fig. 1. Simulated and theoretical CDFs of the regulated test statistics under 
S : Real- valued case, MED with N = 1000. 



B. CND 



(24) For CND, T( C )(X) = h-. Since 71 and j K are asymptoti- 



cally independent, and when N — > 00, ^ 1, 
similar to Theorem 2, we have the following theorem. 
Theorem 5: Under Si, for a fixed K and when N — > 00, 



1. 



■P. (l) 

> C K ' 



where a^- 1 = and the PDF of cjj is given by 



(c) 



(•r) 



9i,i{x + y)gq r , qr (y)dy 



(27) 



(28) 



Now let us look at the case when gi = q r = 1. Since both 71 
and 7ft follow the same Gaussian distribution asymptotically, 
we can easily derive that cj£ ~ AT(0, 4) for real-valued case, 



and c 



(i) 

K 



Af(0, 2) for complex-valued case. 
VI. Computer Simulations 



Computer simulations are presented in this section to val- 
idate the effectiveness of the results obtained in this paper. 
We choose K = 2, and there is one primary signal under Si. 
10000 Monte Carlo runs are carried out in order to compute 
CDF of the test statistics. 

A. Scenario Sq 

We first compare the CDFs of the regulated test statistics 
\/N (X.) — 1) for MED derived using computer simula- 
tions, and theoretic analysis under fixed K assumption of this 
paper and large K assumption of ifTTI . The results for real- 
valued case with N = 1000 are shown in Fig. Q] It is seen 
that the simulated CDF and theoretical CDF derived in this 
paper are close to each other, while the result derived under 
large K assumption are far from the simulated one. 

We next evaluate the accuracy of threshold setting for the 
CND detection using the formula in Section IV. Fig. [2] shows 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Probability of false alarm 



Fig. 2. Comparison of threshold values under So: Real- valued case, CND 
with N = 10000. 



the threshold values at different probability of false alarms 
for real-valued case with N = 10000. For comparison, the 
true thresholds (based on simulations) and the thresholds by 
using the theory for large K in J5] are also included in the 
figure. The proposed theory based on fixed K is much more 
accurate than the theory in and gives the threshold values 
approaching to the true values. 

B. Scenario S\ 

We then compare the detection probability results predicted 
by the formulas derived in this paper and reference (51 with 
large K assumption. The probabilities of detection for CND by 
using simulations and different theoretical formulas are shown 
in Fig. H with N = 10000 and SNR = -15dB (real-valued 
case). For a target probability of false alarm, we first use 
simulations to determine the decision threshold, then obtain 
the probability of detection based on simulations or the related 
formulas. Again, the proposed formula based on fixed K is 
much more accurate than the formula in (5j and gives the 
values approaching to the true values. From Fig. [3] it is also 
seen that, interestingly, for a target probability of false alarm, 
the probability of detection predicted by our method tends to 
be conservative, which seems to be good to primary users in 
terms of protection requirement. Finally, as the results derived 
in this paper are asymptotic, they will approach the simulated 
ones more accurately when the sample size increases. 

VII. Conclusions 

In this paper, theoretic distributions of the test statistics for 
some eigenvalue based detections have been derived for any 
fixed K but large N, from which the probability of detection 
and probability of false alarm have been obtained. The results 
are useful in practice to determine the detection thresholds and 
predict the performances of the sensing algorithms. Extensive 
simulations have shown that theoretic results have higher 
accuracy than existing stuties. 




Probability of false alarm 

Fig. 3. Comparison of estimations for probability of detection under Si: 
Real-valued case, CND with N = 10000 and SNR = -15dB. 
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